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Abstract 

The  Fisher  distribution  is  the  analogue  on  the  sphere  of 
the  isotropic  bivariate  normal  distribution  in  the  plane.  The 
purpose  of  this  paper  is  to  propose  and  analyze  a  spherical 
analogue  of  the  general  bivariate  norral  distribution.  Estim¬ 
ation,  hypothesis  testing  and  confidence  regions  are  also 
discussed. 
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Introduction. 


The  Fisher  distribution  on  the  sphere  is  the  analogue  of 
the  isotropic  bivariate  normal,  that  is  the  Fisher  distribution 
has  circular  contours  of  constant  probability.  However,  in  some 
problems  it  is  desirable  to  have  a  more  general  distribution  on 
the  sphere  with  oval  contours  in  order  to  provide  an  analogue  of 
the  general  bivariate  normal  distribution.  The  purpose  of  this 
paper  is  to  construct  a  suitable  spherical  analogue  of  the  general 
bivariate  normal  distribution  (denoted  as  the  FE ^  distribution 
below). 

After  setting  up  our  notation  in  Section  2,  we  define  the 
8-paraaeter  Fisher-Binghsm  distribution  (FBg)  in  Section  3.  The 
FBj  distribution  will  appear  as  a  5-parameter  sub-family  of  FBg. 
The  limiting  normal  behaviour  of  FB^  for  large  concentration  and 
other  motivating  properties  of  FB^  are  discussed  in  Sections  4  and 
6. 

Sufficient  statistics  for  the  FB^  distribution  are  described 
ir.  Section  3,  estimation  of  the  parameters  in  Section  8,  and  a 
confidence  region  for  the  mean  direction  in  Section  10.  Several 
hypothesis  tests  of  interest  are  discussed  in  Section  9. 

An  example  to  illustrate  the  use  of  the  use  of  the  FB5 
distribution  is  given  in  Section  11.  Analogues  of  the  Fisher- 
Bingham  distribution  in  other  dimensions  are  briefly  mentioned  in 
Section  12. 

Although  the  primary  emphasis  in  this  paper  is  on  the  FB^ 
distribution,  properties  which  applying  to  other  sub-families  of 
FBg  will  also  be  mentioned  where  relevant. 
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2.  Notation 

3  2  2  2 

Let  •  (xcR  :  x'x  •  Xj  ♦  x^  ♦  ■  1)  denote  the  unit 

sphere  in  R^.  We  car,  write  x  in  polar  coordinate!  (6,$) 
defined  by 

Xj  »  cose,  x^  *  sinS  cot*  ,  x^  ■  sing  sin*  (2.1) 

where  0<6<«,  0<«<2«.  If  dx  denotes  Lebesgue  measure  on 

tJj,  then  in  polar  coordinates  dx  ■  sin  6 dd  d*.  Throughout  this 
paper  we  shall  define  distributions  on  3^  terms  of  densities 
with  respect  to  dx. 

A  useful  way  to  plot  spherical  data  is  given  by  Lambert 's 
equal  area  projection  (see  e.g.  Mardia,  1972,  p.215)  defined  by 

*2  •  pcos*  ,  *3  »  paint  (2.2) 

where  p  »  2  sin  (6/2)  ,  0<p<2. 

For  any  matrix  A(n'p),  let  A’  denote  the  transpose  of 
A,  *(j)(nxl)  the  jth  column  of  A,j»l,...,p,  and 
a^(p»l)  the  ith  row  of  A  (written  as  a  column  vector), 
i  •  l,...,n. 

An  orthogonal  matrix  r(3*3)  of  positive  determinant  depends 
on  3  polar  coordinates.  Let  us  denote  by  r  ■  r(*,n,t)  the 
matrix  defined  by 


[  "  C!(l)*  C0*C  *(2)  4  *inC!(3)*  “#inC  !(2)  4  CO*C  “(3): 

(2.3) 


where 

G  • 


cost* 

sin*  cosn 
sin*  sinn 


-sin* 
cos*  cosn 
cos*  sinn 


0 

-sinn 

cosn 


(2.*) 


and  Os*S*  ,  0?n,C<2r  . 

We  next  define  a  concept  we  shall  need  later  in  the  paper. 
Given  a  non-tero  vector  u(3*l)  and  a  spnetric  matrix  A,  let 
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E  •  E(u ,A)  be  a  (3»3)  orthogonal  matrix  such  that  if  v  •  E'u 

•  •  m  m  m  m  ^ 

and  B  •  E'AE,  then 

v2  "  V3  "  °*  *23  "  °*  (2*5) 

Vl>0  •  b22*b33  * 

Call  the  columns  of  E  the  constrained  eigenvectors  of  (u,a). 
Kote  that  e^j  is  proportional  to  u  whereas  e^j  and  e^ 
diagonal ite  A  "as  such  as  possible"  subject  to  being  constrained 
by  *(})*  Also  note  that  defines  a  vector  (whose  sign 

is  determined  by  (2.6)),  whereas  e^  and  e^  only  define  axes 
(whose  order  is  determined  by  (2.6)).  The  constrained  eigenvectors 
can  also  be  viewed  as  the  eigenvectors  after  projecting  A  onto 
the  subspace  orthogonal  to  u  (see  Rato,  1966,  pp  61-62). 

It  is  also  convenient  to  sustsarixe  (u,a)  in  terns  of  the 
site  and  shape  quantities, 

rl  "  W1  "  M  «d  r2  *  b22‘b33»  (2-7) 

respectively,  which  are  invariant  under  orthogonal  changes  of  the 
coordinate  system. 

For  conputstional  purposes,  the  suitrix  E  is  nost  easily 

obtained  by  the  following  two-step  procedure.  First  choose  an 

orthogonal  matrix  H  to  rotate  u  to  the  north  pole 

•  • 

(1,0,0)'.  (In  the  polar  coordinates  of  (2.1)  with 
H  for  T,  choose  6  and  n  ao  that  **•  here  C  is 

arbitrary,  so  for  simplicity  we  csn  take  C  ■  0.)  Then  set 
v«  S'u,  and  C  ■  B'AH.  Secondly,  choose  a  rotation  K  about 

•  •  m  •  •  mm  m 

the  north  pole  to  diagonal ise  C^,  where 

Sh¬ 
is  the  lower  (2*2)  submatrix  of  C.  (Xn  the  polar  coordinates  of 


*22  eli 

*32  *3j 
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(2.1)  with  K  for  f,  take  ;  •  0,  ii  •  0  and  choose  (  to 
satisfy 

tan  •  2^2 3/ ^c22*c33^  * 

ensuring  that  (2.6)  also  holds.)  Then  E  •  HR. 

Note  that  even  after  the  first  stage  the  site  and  shape 
have  staple  interpretations  in  terns  of  v  and  C.  Letting 

-1  » 

ij  and  ^ denote  the  eigenvalues  of  we  have 


'1* 


r2  •  lrl2 


(2.8) 


3.  The  Fisher-Binghas  distribution. 

Define  a  distribution  on  f.«  by  the  density 
3  2 

f(x)  «  ex?!«f  v'x  *  T  6.  (•,  •  x)  }  t  x'x  •  1. 

«.  *  •  j  B  2  J  UK  **  • 


(3.1) 


We  shall  call  (3.1)  the  Fisher-Binghar  distribution  since  the 
first  factor  is  proportional  to  a  Fisher  density  and  the  second 
to  a  Bingham  density.  The  8  parameters  of  (3.1)  are  k;0,  real* 
valued  £2**3'  *  unit  vector  v»  *nd  *r*  orthogonal  matrix 

1  *  [I(l)f:(2)  *  ?(3)3- 

Ve  shall  also  use  the  name  FBg  for  the  full  family  (3.1). 

j  2  f  1  9 

Note  that  the  constraint  T^x*  ■  *  1  taplies  that  a 

tern  *n  th*  axpenent  of  (3.1)  would  be  redundant  in 

the  specification  of  the  density. 

The  family  of  FBg  distributions  is  closed  under  orthogonal 

transformations.  If  x  is  a  random  vector  from  FBg(<(  Bj.Bj.v.r) 

and  8  is  orthogonal,  then  H'x  comes  from  FB8(k,8  ,£  ,H'v,HT). 

•  •  •  •  »  *  *  *  • 

Rote  that  the  transformation  x(*  H'x  can  be  thought  of  as 

•  •  • 

changing  the  frame  of  reference,  with  the  coordinate  axes  in  the 
new  frsae  given  by  the  columns  of  B. 
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Several  interesting  distributions  appear  as  special  cases 
of  FBg.  Besides  the  uniform,  Fisher  and  Binghaa  distributions 
themselves,  ve  also  have  the  following  families,  all  of  which  are 
also  closed  under  orthogonal  transformations. 

(a)  FBfc  (e,6,,fcyr).  Put  v  ■  y^jj  so  that  the  Fisher  axis 
lines  up  with  one  of  the  Binghaa  axes.  Then  the  number  of 
paraaeters  is  reduced  by  2  to  6  parameters. 

(b)  FB.  (*,£,T).  Put  v  •  Yqj  and  set  *  6  say,  with 

£j0.  This  distribution  is  a  5-parar.eter  sub-fanily  of  FB  and 

t) 

is  proposed  here  as  a  spherical  analogue  of  the  bivariate  noraal 
distribution.  The  justification  for  this  proposal  will  be  given 
in  Section  6. 

(c)  FB4  <«,£,  Put  v  *  Y(1)  and  set  “  S3  *  »*?• 

Again,  FB^  is  a  sub-family  of  FB^  but  with  very  different  behaviour 
from  FB^.  Note  that  since  we  cannot  distinguish  between  y^)  *ntJ 
Yqj  here,  FB^  is  only  a  4-paraaeter  family.  If  ,  then 

the  wide  of  the  FB&  density  is  a  small  circle  whose  center  lies 
along  the  axis.  This  distribution  was  introduced  and  studied 

by  Bingham  and  Mardia  (1978) . 

All  of  the  above  families  of  distributions  are  closed  under 
arbitrary  rotations  of  the  coordinate  system.  The  inclusion 
relationships  between  them  are  susnarised  in  Figure  1. 


Figure  1 
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The  Fisher-Bingham  distribution  was  first  introduced  in 

Mardia  (1975,  p.352)  on  the  sphere  il  ,p>2,  It  also  forms  one 

P 

of  a  hierarchy  of  distributions  considered  in  Beran  (1979). 

From  Beran's  point  of  view,  the  Fisher  distribution  contains  an 

arbitrary  linear  function  of  x  in  the  exponent  of  the  density, 

and  FBg  contains  an  arbitrary  linear  and  quadratic  function  of 

x.  The  other  distributions  in  the  hierarchy  include  higher  order 

polynomials  in  the  exponent  of  the  density.  An  extension  of  FBg  to  a 

Stiefel  manifold  was  proposed  in  Mardia  and  Khatri  (1977). 

As  noted  in  Mardia  and  Khatri  (1975),  the  F£.  distribution 

o 

can  be  obtained  by  conditioning  a  trivariate  normal  distribution 
with  arbitrary  mean  vector  and  covariance  matrix  (a  9-parameter 
family)  to  lie  on  the  unit  sphere. 

Unfortunately,  statistical  work  with  the  full  FBg  distribution 
has  been  hampered  by  difficulties  in  estimating  and  interpreting 
the  parameters  (but  see  de  Waal,  1979).  However,  as  we  show  in 
this  paper,  these  difficulties  do  not  apply  to  the  FBg  distribution. 

For  FB^  the  parameters  have  important  and  natural  interpretations, 
and  estimation  is  quite  feasible. 

A.  Limiting  behaviour  of  FBg  for  large  concentration. 

When  the  FBg  distribution  is  highly  concentrated  about  a  point, 
it  is  we 11 -approximated  by  a  bivariate  normal  distribution.  This 
result  generalises  the  well-known  property  of  the  Fisher  distribution 
(see  e.g.  Mardia,  1972,  p.246),  where  an  isotropic  bivariate  normal 
appears.  Details  about  the  closeness  of  this  approximation  in  the 
Fisher  case  ean  be  found  in  Kent  (1978). 

For  convenience  suppose  that  the  orientation  matrix  r  equals 
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I,  the  identify  matrix.  Then  in  the  polar  coordinates  (2.1), 
the  FB^  density  in  (3.1)  takes  the  fora 

g(6,d)  «  exp{x  cos  6  4  gj  sin^  cos2d  ♦  Bj  sin^  sin2d)  •  (4.1) 

Theorem  4.1.  Let  x  -  FB.  (r,  B. ,6_.  1)  and  let  the  parameters 

-  _  o  12. 

x.Bj.Ej  vary  in  such  a  way  that 

*-*•  ,  B 2 "*  d2*  *3^“  **  d3  with  '^<d3sd2<^  ’  <4*2> 

Then  as  x-*«  , 


6  cos  C 

ax* 

~X2 

“o" 

t 

rv"  ,j 

€  sin  d 

A 

0 

—  3 

(4.3) 

where  “  denote?  asymptotically  equal  in  distribution  and  Z 
denotes  convergence  in  distribution. 

Proof.  Using  the  Taylor  series  expansions 
2 

cos  e •  1-6  /2  ♦  •••  ,  sine  ■  B *  • • •  (4.4) 

for  6  small  and  using  (4.2)  we  see  that  (4.1)  is  approximately 
proportional  to 

exp{-J<[62  -2d262  cos2d”2dj62sin*d3)  ,  (4.5) 

which  is  the  form  of  the  limiting  density  in  (4.3).  To  make  this 
argument  rigorous,  it  is  merely  necessary  to  show  that 

(a)  the  approximation  (4.5)  is  adequate  for  |8|<k  ^  and 

(b)  the  probability  mass  associated  with  |6|>k  ®  is  negligible. 
The  details  are  straightforward. 

Similarly,  it  is  straightforward  to  show  that  the  two 


expressions  on  the  left-hand  side  of  (4.3)  are  asymptotically  equal 
in  distribution.  0 
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Corollary  4.1.  Let  x  -  F£j(«,£,I).  If  «r  and  6  vary  in 
auch  a  way  that 

,  6/ic-*d  with  Ojd--i  , 

then 


D  K  |  i  0 

1  t 

2  •  * 

j  0 

1 

•c* 

6  cos  _  . 

j  ■  K 
e  sin  $| 

*2 

x3 

• 

(l-2d)-1 

0 

0 

(143d)'1 

_  _1 

l  *  _ 

l 

— 

5.  Sufficient  statistics 

The  FBg  distribution  forms  a  canonical  exponential  family  in 
its  8  parameters.  If  we  write  a  •  <\>  and  ~  m  T  diag(0, 52.53)^' 
then 

f  (x)  «  exp{o’x  ♦  x*~x}  (5.1) 

and  a  possible  choice  for  the  natural  parameter  vector  is 


(ara2,03’  ^12*  *13  *  T22*  *23*  r33> ' 


(5.2) 


Given  an  (nx3)  data  matrix  X  from  the  FBg  distribution, 

the  corresponding  sufficient  statistic  is 
n 


!  "  n  Ji(xii»xi2**i3*  2xilXi2’  2xilXi3’  Xi2’  2xi2*i3’  Xi3) ' 


(5.3) 

Note  that  t  holds  the  information  contained  in  the  sample  mean  vector 
and  the  sample  dispersion  matrix  about  0, 


x 


—  Z  x.x! 
n  .i.x 


(5.4) 


Unfortunately,  the  other  families  described  above  (FB^.FB^  and  FBg) 
are  not  canonical  exponential  families,  but  instead  from  curved 
exponential  families.  In  each  case  the  parameter  space  has  fewer 
than  8  distensions,  but  the  minimal  sufficient  statistic  is  still 
given  by  (5.3). 
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These  remarks  about  FB^.FB,.  and  FB^  apply  only  when  all 
the  parameters  are  unknown.  A  neater  situation  arises  if  we 
suppose  that  the  mean  direction  v  (and  also  possibly  the 
concentration  k)  is  known.  With  this  knowledge  FE^,  FB^  and 
FB^  (and  also  FBg)  now  become  canonical  exponential  families. 

For  definiteness  we  consider  FB,..  After  rotating  the 
coordinate  system  so  that  v  «  becomes  equal  to  (1,0,0)*  the 

density  is  proportional  to 


2  2 

f  (x)  *  exp{‘ x,  4i*:<2~X3>  *  2i2  *23^  ’ 


(5.5) 


or  in  polar  coordinates. 


g(P,$)  “  exp{*  cos  6  *  t  sin  6  cos2fe—A))  , 


(5.6) 


where 


• £  cos  2X  ,  f2  ■ £  sin  2X  . 


(5.7) 


The  parameter  X€[0,2t)  describes  the  direction  of  the  major 
axis  -  see  Section  6.  Then  the  natural  parameter  and  sufficient 


statistic  are  given  by 


(■'•W*  *nd  n  i^1C*il*xi2"  *i3’  2xi2Xi3)f  *  (5,8) 
If  k  is  also  known,  the  natural  parameter  and  sufficient 


statistic  become  slightly  simpler;  namely 

<W  -nd  n  1j1(xi2_xi3*2xi2xi3>  ’  * 


(5.9) 


6.  Properties  of  the  FB^  distribution. 

The  FB^  density  was  defined  in  Section  3  by 

f(x)  «  exp UY(2£  *  ®^Y(2)X^2"^(3)^2^  * 

As  we  shall  see  below  the  parameters  can  be  described  as 
follows:  r* 0  is  the  concentration,  0*0  describe  the  ovalness. 


10 


Y(l)  *S  106130  direction  or  pole,  ’(-7)  is  the  niajor  axis,  and 
^(3)  *s  t*1e  m*nor  axis.  Note  that  y^  and  y^  are  determined 
only  up  to  sign,  bo  that  they  do  indeed  define  axes  rather  than 
directions . 

If  we  rotate  to  the  frame  of  reference  defined  bv  the  columns 
of  T,  the  density  f(x)  takes  a  particularly  simple  form.  For 
this  reason  we  shall  call  this  transformation,  x(*  x*  *  r'x, 
the  transformation  to  the  population  standard  frame  of  reference.  The 
density  for  x*  takes  the  form 

i(x*)  =  exp{^  x*  ♦  6(x*2  -  x*-)} 
or  in  polar  coordinates 

2 

g(S,$)  «  exp{<  cosS  +  £  sin  e  cos2p}  . 


A  sample  analogue  to  the  population  standard  frame  will  be 
defined  in  Section  8. 

As  stated  in  the  introduction,  the  FB^  distribution  is 
proposed  here  as  a  spherical  analogue  of  the  bivariate  normal 
distribution.  The  following  properties  show  why  this  is  a  sensible 
proposal.  We  need  to  suppose  here  that  2t<<  to  ensure  the  correct 
behaviour. 


(a) 

(b) 


FB^  and  the  bivariate  normal  are  both  5-parameter  families. 
The  contours  of  constant  probability  near  the  pole 
are  approximately  ellipses  with  major  and  minor  axes  y^) 
and  respectively.  (This  property  follows  easily 

from  (4.4).) 


(c)  The  geometric  average  of  g(6,$)  over  circles  of  constant 
latitude  is  proportional  to  a  Fisher  density,  that  it. 
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f2” 

log  g(6,t)<H  •  « cose  ♦  constant. 

Jo 

Thus,  in  this  sense  FB^  is  a  natural  extension  of  the  Fisher 
distribution. 

(d)  As  6  goes  .from  0  to  *  for  fixed  6,  g(6,4)  decreases 

monotonical ly.  Thus  g(E,<)  is  unioodal  on  all  great  circles 
through  the  pole. 

(e)  For  large  values  of  the  concentration  parameter  k,  FB,.  is 

approximately  the  same  as  a  bivariate  normal  distribution 
with  mean  and  major  and  minor  axes  *nt*  Y(3) 

respectively.  (See  Corollary  4.1.) 

Note  that  the  larger  FB^  family  is  not  a  suitable  spherical 
analogue  of  the  general  bivariate  normal  distribution  because  it 
has  one  too  many  parameters.  In  Theorem  4.1  we  saw  that  this 
"extrrf' parameter  is  "asymptotically  unidentifiable"  for  large 
concentration. 

Hence  ve  have  introduced  a  further  constraint  (S^* - 6^)  to 
define  the  FB^  family.  To  some  extent  this  constraint  is  arbitrary 
(in  fact  a  different  constraint  was  proposed  in  Kent,  1980).  However, 
the  constraint  used  here  does  have  some  theoretically  attractive 
properties  ((c)  and  (d)  above).  Moreover,  as  we  shall  see  in  the 
next  section,  with  this  definition  of  FB^  the  normalization  constant 
takes  a  reasonably  tractable  form. 

7.  Moments  of  FB, . 

So  far  ve  have  not  dealt  with  the  normalization  constant  of 
the  FB^  distribution, 

e(*,6)  *  f  f  *XP{* eotS  *  8  «i*»2e  eos2$}sinB  d$  de.  (7.1) 
Jo  Jo 
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Using  the  results 

»ii  /2 

I  (sin?)a(cos£)b  de  -  iB(^i  ,  ,  (7.2) 

-o  z  1 

(Abramowit2  and  Stegun,  1972,  (6.2.1),  p.256)  where  B(*,*>  is 
the  beta  function,  and 

I  e*0080  sin2ved6  -  air(v*J)(J<rv  I  (<)  (7.3) 

'o  v 

(Abramowitz  and  Stegun,  1972,  (9.6.16),  p.376)  where  1  (*c)  is 

V 

tne  modified  Bessel  function,  we  can  expand  c(k,£)  in  a  scries 


c(<,£) 


t  rJk+i l 

L  7(k*l)  C  }  I2k-i 


(7.4) 


Since  sequences  of  Bessel  functions  can  be  quickly  and  easily 
computed  by  the  method  of  Amos  (1974),  formula  (7.4)  provides  a 
quick  and  simple  method  of  calculating  c  (*,£).  Note  that 
c(o,o)  ■  4r,  the  surface  area  of  the  sphere,  and 
c(*,o)  •  4t<  *  sinh k  , 

the  normalizing  constant  for  the  Fisher  distribution. 

For  large  k  (with  26/*<l  fixed)  we  have  from  Corollary 
4.1  the  asymptotic  formula 

c(k,£)  •  2reK[(K-26)(^2c)]"1/2  .  (7.5) 


Consider  a  random  vector  x  -  FB^(k,£,I)  , 
(7.1)  under  the  integral  sign  and  writing  c  • 
etc.,  we  find 

2 

Ex.  •  c  /c  ,  Ex,  ■  c  /c 
Ik  1  ric 

E(x2-x2)  •  Cg/c  , 

2  2  2 

and  hence  since  Xj+x^+Xj  •  1, 


Differentiating 
c(k,6),  ck  »  8c(k,B)/ci<, 


(7.6) 


SB!  gggggggfgg*?1 


13 


•> 

Lx*  ■  (c-c  ♦  c ,)/2c 

2  “  6  V.r, 

Ex^«  (c-c  -  c  J/2c. 

J  Kk  f 

For  later  use  write 

n«  E*J  and  o?*  Ex?,  j  -  1,2,3.  (7.8) 

By  symmetry,  most  of  the  other  moments  of  interest  equal  0, 
Kx^)  «  E(Xj)  ■  0,  (7.9) 

E^X2*3>  *  °»  (7.10) 

E(x2x2)  »  E(XjX^)  •  0.  (7.11) 


8.  Moment  Estimation  for  FB  . 

- -  3 

Let  x  ...... x  be  a  sample  from  FB.(*,£,T).  The  standard 

« i  *r»  d  % 

way  to  estimate  the  parameters  of  FB^  is  to  use  maximum  likelihood 
•  * 

estimates  r,£,F.  However,  it  does  not  seem  possible  to  obtain 
explicit  expressions  for  the  m.l.e.s,  so  iterative  methods  must 
be  used  to  find  them. 

In  this  section  ve  propose  simpler  estimates  which  we  call 
the  moment  estimates  «,£,r  for  the  parameters  of  FB..  They  have 
the  following  properties. 

(a)  The  moment  estimates  are  consistent  estimates  of  the 

true  parameters  and  hence  provide  suitable  starting  values 
for  maximum  likelihood  iteration. 

(b)  The  orientation  matrix  T  can  be  calculated  explicitly. 

(c)  If  either  the  eccentricity  2£/ic  is  small  (the  usual  case 
in  practice)  or  if  x  is  large,  then  the  moment  estimates 
are  close  to  the  m.l.e.s. 

(d)  If  the  data  is  highly  concentrated,  the  concentration 

»  m 

parameters  x,g  can  also  he  calculated  explicitly. 
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More  specifically  the  moment  estimates  axe  defined  as 
follows.  Let  x  and  S  be  the  sample  mean  vector  and  the 
sample  dispersion  matrix  about  0  as  in  (5.4).  Then,  in  the 
terminology  of  Section  2,  T  is  defined  to  be  the  matrix  of 
constrained  eigenvectors  for  (x,S).  Further,  letting 
r^  ■  jjxjj  and  r^  denote  the  size  and  shape  quantities  for 
(x,S),  the  estimates  k,£  of  concentration  parameters  are 
determined  implicitly  by  the  equations 

rj-c^/c  ■  0,  r2“ce/c  ■  0.  (6.1) 

For  large  <,  the  use  of  (7.5)  leads  to  the  explicit  solution 


<  -  (2-2r1-r2)'1  ♦  (2-2r1*r2)‘1 
6  -  i((2-2r1-r2)”1  -  (2-2r^r2)"1 }  . 


(6.2) 


The  orientation  matrix  T  has  been  chosen  so  that  for  x* 
the  sample  analogues  of  (7.9)-(7.10)  (but  not  (7.11))vill  hold. 

For  this  reason  ve  shall  term  the  transformation  xk  x*  ■  T ' x 
the  transformation  to  the  sample  standard  frame  of  reference. 

The  rationale  behind  the  moment  estimate  T  is  as  follows. 

The  first  column  is  the  mean  direction  of  the  sample,  which 

is  also  the  m.l.e.  of  the  mean  direction  under  a  Fisher  distribution. 

If  the  eccentricity  28/k  is  not  too  large  (which  is  the  most 
important  case  in  practice),  then  y^jj  will  also  be  close  to  the 
mJ.e.of  the  mean  direction  for  FBj.  Further,  if  the  true  mean  direction 
Yqj»y^jj  were  known,  then  (7.3)-(7.4)  would  ensure  that  y^ 
and  Y(3)  would  be  the  m.l.e.s  of  the  major  and  minor  axes,  respectively. 
Bence  if  26/k  is  not  too  large,  the  moment  estimates  should  be 
nearly  as  efficient  as  the  maximum  likelihood  estimates. 


A  similar  situation  arises  for  large  «  tritan  FB^  is  close 
to  a  bivariate  normal.  Then  the  moment  estimates  and  m.l.e.s  for 
FB^  will  both  be  close  to  the  corresponding  m.l.e.s  for  the 
bivariate  normal. 


9.  Some  hypothesis  tests. 

Zn  this  section  we  describe  several  large-sample  hypothesis 
tests  of  interest.  All  of  these  tests  can  be  carried  out  using 
the  following  general  result. 

Theorcr  9.1.  Consider  n  independent  identically  distributed 

observations  from  a  model  H.  with  parameters  (r,>.)  of  dimensions 

•  •  • 

p  and  q  respectively,  and  consider  a  null  hypothesis  Ho  :  >.  •  0. 
Suppose  that  under  H^,  the  model  forms  a  canonical  exponential 
family  for  t  with  minimal  sufficient  statistic  u,  and  that  under 
Hj  (with  -  known)  the  model  fores  a  canonical  exponential  family 
for  }.  with  minimal  sufficient  statistic  w.  Define  Kao's  score 

m  m 

statistic  by 

V  •  [w-E(wJu)]'  Var(w|u)~*  [w-E(w|u)D  t  (9.1) 


*  A 

where  ell  moments  are  calculated  under  H  with  *•*  («  being  the 

o  .  . 

m.l.c.  of  t  under  H  ). 

o 

Then  asymptotically  as  the  sample  aise 


and  further  Vu  is  asymptotically  equivalent  to  the  likelihood 
ratio  statistic  -2  log  l  for  I  vs  I. .  We  rejeet  I  if  V  is  too  large. 

P  A  •  • 

Fr oof.  This  result  is  a  special  case  of  a  general  result  in  Cos  ond 
linkley  (1974),  p.324,  equation  after  (S.i).  lee  also  Kao  (1973),  p.418 

D 

Rote  that  the  score  statistic  WL  is  nsually  simpler  to 

m 

calculate  than  the  likelihood  ratio  statistic  because  only  a 
parametric  model  under  i#  need  be  fitted.  Seme  hypothesis  tests 
which  fit  into  this  framework  mill  mow 


be  described 
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(a)  H  :  Fithcr  vs  H.  :  FB. 

O  A  J 

A  Fisher  distribution  with  concentration  parameter  k  and 
mean  direction  v  forms  a  canonical  exponential  family  vith 
natural  parameters  («Vj,rvj,«Vj)  and  sufficient  statistic 

u  ■  n~*£x..  The  m.l.e.s  satisfy 

v  •  n"1  J x^ /r j  •  ■  tj  (9.3) 

where  r.  •  Hn”*  Jx.  ||  is  the  resultant  length,  also  used  in  (8.1). 

Kou  let  I  be  an  orthogonal  matrix  which  depends  on  the  data  only 

* 

through  v,  and  whose  first  column  is  given  by  h^  •  v.  Let 
y.  •  H'x.  ,  i*l,...,n.  Then  from  (5.9),vith  *  •  e  and 

•  t  .  .. 

•  v  assumed  known,  the  model  under  Hj  fores  a  canonical  exponential 
family  vith  sufficient  statistics 

vi  *  £i<»u**i»>  •  w2  •  loir’ll  ■ 

Nov  the  assumption  that  the  x^  come  from  a  Fisher  distribution 

F(k,v)  is  equivalent  to  the  assumption  that  the  y.  come  from 
•  •>* 

F(«, (1,0,0)*).  Further,  for  a  random  vector  y  *  F(i, (1,0,0) *),  we 
have  by  symmetry 

*(y2“7§>  •  0  .  I(y2y3>  -  0 

*(fJf,>  •  *<f2fJ>  •  0 

KFjfJ-FjfJ)  •  0  ,  *<FjF2Fj)  •  0  ,  J- 1,2.3, 
and  me  have  by  (7.2)  and  (7.2), 

Kyj-yJ)  •  *(*yjyj)  •  (i/2 >“2  («>  . 

■sees  the  test  statistic  (t.l)  takes  the  fen 

wV1  ,  ,  , 

*  tyl)  -  (*i  *  ”i>  -  *a  • 


(t.3) 
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Using  Che  notation  of  Section  2,  it  is  assy  to  check  that 

2  2  ? 

Vj  ♦  equals  the  squared  shape  quantity  r^  for  (x,S).  In 

particular  it  is  clear  that  the  statistic  W’u  does  not  depend 

on  the  arbitrariness  in  the  choice  of  the  second  end  third  columns 

of  H. 

In  practice  this  test  Bight  be  used  in  the  following  situation. 

Given  a  set  of  spherical  data,  an  experimenter  Bight  first  look 

for  directionality  by  testing  H  :  uniform  vs  H, :  Fisher  (the 

o  i 

Rayleigh  test).  If  this  null  hypothesis  is  rejected  he  might 
assess  the  circular  syraetry  of  the  data  about  the  pole  by  using 
the  test  described  here. 

For  large  concentration,  this  test  reduces  to  a  test  of 
sphericity  for  the  bivariate  normal  distribution;  see  for  example 
Mardia,  Kent,  and  Bibby  (1979)  p.  134. 

(b)  Hq:  Fisher  vs  FBg. 

As  in  (a),  this  is  a  goodness-of-fit  test  for  the  Fisher 
distribution,  but  here  with  a  more  general  alternative  in  mind. 

The  calculations  are  sinilar  to  those  in  (a)  but  sonewhat  more 
involved.  Details  are  given  in  Mardia  and  Holsws  (1980). 

The  analogous  test  on  the  circle  was  given  by  Cox  (1975);  see 
also  Cox  and  Barndorff-Nielaen  (1979),  p.291. 

(c)  *ot  Bingham  vs  Bji  FB#. 

This  test  is  included  here  beceuse  it  fits  the  assumptions 
of  Theoren  9.1,  hut  the  application  of  this  test  is  soaevhat 
different  froa  (a)  and  <b).  Suppose  an  experimenter  hss  data  to 
uhieh  he  would  like  to  fit  a  Binghaa  distribution.  The  data  are 
suspected  to  be  hot  not  known  to  be  antipodally  symetric.  There 
are  two  ways  to  proceed  in  this  situation. 

<i>  First  use  a  Kayleigh  ^test  Q*  i»«  *  X*>  to  test  unifornjM^ 


IK* 
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Fisher.  If  uniforsity  is  accepted,  then  antipodal  sytmetry  can  be 
presumed  and  a  Bingham  distribution  can  be  fitted. 

(ii)  First  fit  a  Bingham  distribution,  and  then  assess  the  goodness- 
of-fit  by  using  the  test  of  this  section.  This  latter  approach  seems 
more  suitable  when  the  data  is  known  at  the  outset  not  to  be 
uniformly  distributed. 

Under  H^,  5,  the  sample  dispersion  matrix  about  0,  is 
sufficient  for  the  parameters.  If  x  comes  from  a  Bingham 
distribution,  we  have  by  symmetry 

£(>^1  ■  0  ,  E(x.x  x^  •  0,  i, j,k  ■  1,2,3 

so  that 

E(x)  -  0  ,  E(xiS)  •  0  . 

Further,  since  the  m. l.e.s  of  the  Bingham  parameters  are  chosen 
so  that  E(xx')  -  S,  ve  have 

Var(xjs)  ■  Var(x)  •  ^  S.  Hence  the  score  statistic  takes 
_  _  .  n  . 

the  form 

n  x’S^x  -  X3  •  (9.6) 

•*  Kate  that  when  the  data  is  uniformly  distributed,  S  s  »I,  so 
that  (9.6)  reduces  to  the  Rayleigh  test  statistic. 

10.  A  confidence  interval  for  the  mean  direction. 

Let  *!»•••.*,  be  ■  sample  from  FB5(x,8,D.  First,  transform 

*0  the  population  standard  frame  (see  Section  6)  y.  -  r*x.,iwl . n, 

•  *1 

eo  that  the  yj  form  a  sample  from  FB^c.B.I).  Then  by  the  central 

limit  theorem  and  (7.6)-<7.ll)  we  aee  that  the  sample  mean 

•  «  lI?i  **  asymptotically  normally  distributed  with  mean  vector 

—1  2  2  2 

(p,o,o)*  and  covariance  matrix  n  diag(«|t0j,c3). 


A 
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The  sample  seen  direction  for  the  y.  (which  is  also  the 

moment  estimate  of  the  true  seen  direction  of  the  y.)  is  defined 

& 

m 

by  •  y / J’  y ||  .  Consider  the  two  coordinates  (*21,Y31*'  °* 

Than  by  a  general  result  on  transformations  (see,  e.g.  Rao, 

m  m 

1973,  p.  3S7),  (■y^j.Yjj)  is  also  asymptotically  normal ly 
distributed. 


Hence,  an  asymptotic  100(l-u)J  probability  region  about  the 
scan  direction  of  F£^(«,f,r)  is  given  in  population  standard 
coordinates  by  the  ellipse-like  region 

(Y(1jCf«3  :  Yj^O,  n  w2(y2^=2  *  Y31/c3>*  x2;o}  *  (10.1) 

2  2 

where  x,.  denotes  the  upper  o  critical  value  of  the  x 
*5®  2 

distribution. 

Inverting  (10.1)  provides  a  confidence  region  for  the  true 
scan  direction  about  the  saople  seen  direction.  More  specifically, 
suppose  x.,...,x  cose  from  FB.(k,B,D  with  true  mean  direction 

w  A  w  • 

v  -  Y(i)>  Let  r  be  the  moaent  estimate  of  T  (sec  Section  7)  and 
transform  to  the  sample  standard  frame,  x*^  •  T*  x^,  i*l,...,n, 
and  let  v*  ■  r*v.  Then  an  asymptotic  confidence  region  for  v* 
is  defined  by  the  ellipse-like  region  on  the  sphere 

<v*eOj  *  v*>0,  n  u2(vJ2/o2  ♦  vJ2/oJ)<  x*.,)  (10.2) 

2  2 

Of  course,  in  practice  end  Oj  will  mot  be  known,  but  must 

be  replaced  by  any  consistent  estimates,  (hie  possibility  is  to 
estimate  «  end  0  end  then  use  (7.S)-(7.7).  lewever,  a  simpler 
technique  is  to  just  use  the  sample  moments 
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u  •  n  •  Cj  •  n  *1**2  *nd  °3  "  I  **3  •  (10*3> 

respectively. 

Note  that  in  moving  fro*  (10.1)  to  (10.2)  we  have  switched 
from  s  frame  of  reference  about  the  true  *ean  direction  to  a  frame 
about  the  sample  scan  direction.  However,  since  these  two  points 
lie  within  0(n  *)  of  one  another  in  probability,  the  complications 
arising  from  this  switch  are  negligible  and  (10.2)  remains 
asymptotically  valid. 

When  using  the  equal-area  projection  (2.2)  (in  the  sample 
standard  frame)  and  when  using  the  estimates  (10.3),  a  confidence 
region  asymptotically  equivalent  to  (10.2)  is  given  by 

{(*2.*3)  :  n  ;2(x2/Oj  ♦  tJ/Cj)  <  x*;o).  (10.4) 

11.  Example. 

Creer,  Irving  and  Ksirn  (1959)  measured  directions  of  magnetism 
at  n *  34  sites  in  the  Great  Whin  Sill.  Their  data  is  summarized 
in  Table  1,  column  (b)  of  that  paper,  pp. 311-312  (excluding  sites 
32  and  34).  Their  use  of  declination  (D)  and  inclination  (1)  is 
related  to  our  use  of  polar  coordinates  in  (2.1)  by  6  ■  90  *1  , 
i  •  360°  -D.  [Change  program  aecordinly). 

The  summary  statistics  ere  given  by 


'  0.083“ 

'0.045 

-0.075 

0.014” 

-0.959 

*  • 

m 

■0.075 

0.921 

-0.122 

-0.131 

0.014 

-0.122 

0.034 

_  — 

- 

- 

fro*  which  we  find  r^  •  0*971  ,  tj  •  0*0229. 

The  *o*ent  and  maximum  likelihood  estimates  of  location  are 
both  given  to  three  decimal  places  by 
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0.165 


r 


i  T 


0.085  -0.976 

-0.987  -0.108  -0.117 

0.134  -0.172  -0.976 


and  the  estimates  of  concentration  by 

*  •  42.16  6  •  9.27  (exact  moment) 

k  *  41.76  6  •  8.37  (asymptotic  moment) 

*  *  42.16  6  ■  9.28  (maximum  likelihood). 

The  data  and  a  952  confidence  region  for  the  mean  direction 
baaed  on  (10.2)  and  (10.3)  are  given  in  Figure  2,  plotted  using 
the  equal-area  projection  of  (2.2)  in  the  sample  standard  fraee. 

The  hypothesis  test  of  Section  9(b)  and  the  corresponding 
likelihood  ratio  test  yield  the  values 

«u  -  5.96  and  -21ogL  •  6.55. 

Since  the  upper  52  critical  value  of  x2  5.99,  both 
statistics  show  moderate  evidence  of  a  departure  from  a  Fisher 
distribution. 

In  this  example  an  alternative  approach  was  used  by  Creer  et.al. 
They  managed  to  transform  the  data  to  approximate  circular  symmetry 
and  then  to  use  statistical  techniques  applicable  to  the  Fisher 
distribution. 


12.  Analogues  in  other  dimensions. 

Much  of  the  theory  in  this  paper  extends,  at  least  in  principle, 
to  other  dimensions.  The  analogue  of  the  full  Fisher  *ingham 
family  on  the  unit  sphere  in  R*,  p>2,  can  be  written  in  the 
general  form  (2.1)  with  the  sunset ion  from  2  to  3  replaced  by  a 
summation  from  2  to  p  (as  in  Reran,  1979). 

An  analogue  of  Fl^  can  be  obtained  by  introducing  the  constraint 
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On  the  circle,  p*2,  this  analogue  of  FB^  is  no  more  general  than 
the  von  Mises  distribution  itself.  For  general  p>3,  provided 
|6j j  <  </2,  j  ■  2, . . .,p,  properties  analogous  to  those  in  Sections 
valid.  When  p>3,  moment  estimation  still  can  be  carried 
although  the  normalization  constant  seems  to  become  more 
complicated  as  the  dimension  increases. 


Figure  2:  Floe  of  the  Greet  Vhin  Sill  data  using  an  equal>area 
projection  in  the  sample  standard  franc t  and  a  952  confidence 
•llipee  for  the  aean  direction. 
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